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INTRODUCTION 


As  performance  requirements  for  rocket  motors  are  advanced,  it 
becomes  increasingly  important  to  understand  and  control  the  influence 
of  defects  in  the  motor  grain.  One  of  the  most  potentially  serious 
types  of  defects  are  cracks  which  may  occur  during  curing  or  later  in 
storage.  These  cracks  may  act  as  thermal  insulators  and  may  signifi¬ 
cantly  alter  subsequent  motor  performance.  Attempts  to  develop  and 
validate  mathematical  models  for  predicting  behavior  of  such  cracks  in 
motor  grain  have  been  hampered  by  lack  of  an  extensive  data  base  and 
lack  of  quanti f ’cation  of  three  dimensional  effects. 

In  a  recent  study  [Ref.  Ij,  some  progress  was  made  towards 
development  of  a  aata  base  for  understanding  the  near  tip  behavior  of 
cracks  in  inert  propellant  samples.  These  studies  revealed  that  a 
significant  amount  of  severe  blunting  developed  at  the  crack  tip  prior 
to  crack  extension  but  that,  in  the  early  stages  of  crack  opening,  a 
reduction  in  the  order  of  the  stress  singularity  occured  when  the  crack 
intersected  the  free  surface  of  the  inert  propellant,  and  this  reduction 
was  accurately  measured  at  the  free  surface  by  a  grid  method  which  did 
not  influence  the  near  tip  behavior.  The  variation  in  the  stress 
singularity  order  with  direction  normal  to  the  free  surface  was  then 
studied  using  four  point  bend  frozen  stress  photoelastic  specimens. 

All  of  the  above  studies  were  conducted  on  "artificial"  or  man-made 
straight  front  cracks  and  seemed  to  indicate  that  the  binder,  though 
only  present  in  small  volume  fraction,  controlled  the  near  tip  process. 

Since  it  is  well  known  that  through  cracks  "thumbnail"  when  they 
grow  in  plates  of  continua  such  as  metals  and  plastics  under  Mode  I 
loading,  and  that  this  changes  the  boundary  intersection  angle  between 
the  crack  front  and  the  free  surface  of  the  specimen,  the  question  arose 
as  to  what  the  effect  of  thumbnailing  would  be  upon  the  free  surface 
singularity  order,  its  distribution  through  the  thickness  and  the  near 
tip  displacement  and  strain  fields.  Also  at  issue  was  the  influence  of 
a  surface  layer  of  binder  on  specimens  of  inert  propellant  tested  and 
reported  on  in  Ref.  1.  In  order  to  clarify  these  issues,  the  following 
program  was  developed  and  conducted. 

PROGRAM  OBJECTIVES 

The  objectives  of  the  present  program  were: 

To  measure  the  components  of  the  displacement  field  in  the  free 
surface  of  specimens,  at  intervals  of  time  during  constant  head  rate 
separations  producing  extension  of  the  crack  in  "biaxial"  specimens  of 
inert  propellant  without  surface  layers  of  binder  and  specimens  of  pure 
binder  and  extract  surface  strain  distributions  from  the  data. 

To  estimate  the  order  of  the  stress  singularity  or  dominant 
eigenvalue  for  the  crack  at  the  free  surface  for  a  growing  crack  and  to 
study  the  crack  front  shape  and  its  influence  thereon. 
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PROGRAM  SCOPE 


The  program  consisted  of  three  phases:  1)  Fabrication  of  cracked 
biaxial  specimens  of  both  inert  propellant  and  pure  binder  (with 
additives),  2)  Measurement  of  displacement  fields  near  crack  tips  and 
determination  of  near  tip  strain  fields  and  singularity  order  at  the 
free  surface  for  grown  cracks  along  with  measurement  of  crack  front 
shapes  for  both  inert  propellant  and  pure  binder.  Determine  singularity 
order  distribution  through  thickness  using  the  frozen  stress  method  on  a 
four  point  bend  specimen  for  naturally  grown  "thumbnai led"  cracks. 

3)  Preparation  of  a  final  comprehensive  technical  report  describing  the 
results  of  the  effort. 

WORK  PLAN  AND  ACCOMPLISHMENTS 

The  program  conducted  can  be  conveniently  separated  into  two  parts: 

I)  Mode  I  constant  head  rate  tests  on  "biaxial"  cracked  specimens 
of  pure  binder  and  inert  propellant  with  measurements  of  surface  dis¬ 
placements  and  crack  front  shapes  during  growth  followed  by  extraction 
of  near  tip  strain  fields  and  the  free  surface  stress  singularity  or 
dominant  eigenvalue  order,  and  2)  Frozen  stress  tests  on  four  point  bend 
specimens  containing  thumbnai led  natural  crack  fronts  to  measure 
variation  of  the  lowest  dominant  eigenvalue  with  direction  normal  to  the 
free  surface. 

BIAXIAL  TESTS;  METHODS,  MEASUREMENTS  AND  RESULTS 

Preliminary  tensile  tests  on  the  pure  binder  yielded  a  value  of 
Elastic  Modulus  E  =  112.3  psi  =  0.776  MPa  and  Poisson's  Ratio  v  =  0,47 
at  room  temperature  below  a  nominal  global  strain  (ey)  of  *  20%. 

The  "biaxial"  test  specimen  configuration  is  pictured  in  Figure  I 
and  is  the  same  as  employed  in  Ref.  1.  The  pure  binder  material  con¬ 
sisted  of  polybutadiene  rubber  with  chemical  additives  used  to  enhance 
its  capability  to  bond  to  propellant  particles.  It  was  manufactured  by 
Morton  Thiokol  under  the  name  R45M  Gumstock.  The  inert  propellant  con¬ 
sisted  of  this  binder  with  embedded  polyhedral  particles  of  potassium 
chloride,  the  latter  occupying  a  volume  fraction  of  approximately  75%. 
The  particle  sizes  ranged  from  .01  to  0.5  mm.  The  inert  propellant  was 
also  supplied  by  Merton  Thiokol  under  the  name  H-24  inert  propellant. 

A  grating  consisting  of  squares  0.2  mm  on  a  side  was  applied  to  the 
inert  propellant  specimen  using  the  following  procedure: 

•  Spread  titanium  dioxide  powder  over  the  specimen  surface  near 
the  crack  tip. 

•  Apply  a  mixture  of  silicone  grease  and  titanium  dioxide  to  the 
surface. 

•  Place  a  metal  screen  in  the  surface  mixture  and  remove  excess 
mixture. 


2 


1 


DIMENSIONS  IN  MILLIMETERS:  Initial  Crack  Length  (2a)  ^38  mm 

Specimen  thickness  =  15.25  mrn 

^ig.  1  Biaxial  Test  Specimen  Dimensions 

•  Evaporate  aluminum  over  the  screened  surface. 

•  Remove  the  metal  screen. 

The  same  size  grating  was  produced  on  the  binder  using  the  follow¬ 
ing  procedure: 

•  Apply  a  special  adhesive  to  the  binder  surface. 

•  Place  a  metal  screen  on  it. 

•  Evaporate  aluminum  on  this  surface. 

•  Remove  the  metal  screen  as  before. 

The  total  thickness  of  the  resulting  grating  was  less  than  2.5x10"^ 
mm. 


Tests  were  conducted  at  constant  head  rates  of  2.5  mm/min  and  at 
25.4  mm/min  and  photographs  of  the  grid  region  were  taken  at  various 
time  intervals.  A  schematic  of  the  photographic  setup  is  shown  in 
Figure  2.  The  lighting  provided  by  the  ring  illuminator  and  reflector 
was  critical.  Figure  3  shows  typical  load-extension-time  records  for 
tests  on  both  propellant  and  pure  binder  at  each  rate  showing  both  the 
strengthening  and  stiffening  effect  of  the  rigid  particles.  Figure  4 
illustrates  a  typical  set  of  raw  data  for  an  inert  propellant  test  and  a 
pure  binder  test  illustrating  the  processes  involved  in  crack  extension. 
At  a  given  instant  during  crack  growth,  a  dark  powder  was  "shot"  into 
the  crack  opening  to  mark  the  crack  front  shape.  A  summary  of  the 
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SPECIMEN 


Fig.  2  Schematic  of  Photographic  Setup 


biaxial  tests  conducted  and  the  corresponding  head  rates  and  dominant 
eigenvalue  data  zones  are  given  in  Table  1. 


Data  from  photographs  such  as  those  in  Figure  4  were  digitized  and 
displacement  and  strain  fields  were  extracted  following  the  procedures 
of  Ref.  1  (Appendix  A).  Typical  results  are  shown  in  Figure  5.  Examin¬ 
ation  of  results  revealed  that  the  effect  of  changing  the  head  rate  by 
an  order  of  magnitude  altered  the  displacement  fields  (and  also  the 
strain  fields)  but  the  distribution  contours  were  of  the  same  general 
form.  Consequently,  we  shall  present  digitized  data  here  only  for  the 
slow  rate.  Figure  5  presents  a  complete  set  of  digitized  results,  (U, 

V,  e  ,  e  ,  y  )  for  a  global  strain  level  of  e  =  4%  from  an  inert  pro- 
x  y  xy  y 

pellant  (]_P)  test  (Figs.  5a  through  5e)  and  a  set  for  a  global  strain 

level  of  ty  =  4.7 %  from  a  pure  binder  (PB)  test  (Figs.  5f  through  5j) 

for  purposes  of  comparison.  Comparing  the  displacement  fields  in  Figs. 

5a  and  5b  with  those  in  Figs.  5f  and  5g,  respectively,  we  find  similar 

global  contour  shapes.  In  determining  e  and  e  ,  data  were  averaged 

x  y 

from  the  two  sides  of  the  crack  plane  to  eliminate  any  possible  in-plane 
shear  mode  effect. 


Load  (lb) 


Load-Extension-Time  Records  -  Inert  Propellant  (IP)  and  Pure  Binder  (PB 


(c)  e  =  4.76% 

i  y 


Fig.  4 


Crack  Opening  and  Growth  in  Specimen  P-2 
a)  Opening  5)  Blunting  c)  Void  Formation 

d)  Coalesence  of  Void  with  Main  Crack 

e)  Crack  Extension 
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TABLE  1 


ui.1x1.3l  Test  Data  for  »  Determination 

u 

Grating  Frequency  5°. /mm 


Test  No. 

Held  Pate  mm/m  in 

Data  Zone  for 

Inert 

Prope 1 1  ant 

P-1 

2.5 

3. 5-5.0 

P-2 

2.5 

I. 5-5.0 

P-3 

25.4 

1.5-4. 5 

P-4 

25.4 

1. 5-5.0 

Pure 

Binder 

8-1 

2.5 

1. 0-6.0 

B-2 

2.5 

1. 5-5.0 

B-3 

25.4 

I. 0-6.0 

B-4 

25.4 

1. 5-5.0 

B-5 

2.5 

1  5-5.0 

B-6 

2.5 

1. 5-5.0 

Computer  plots  of  tnj 

in  Figs.  5c  througr  5,?  may 
5h  throuy'  5j.  The  icur:  <i 
sketched  in  for  purposes  jf 
strain  levels  greate-  them 
definitions  //ere  used  ner-?. 


strain  distributions  r  ,  -  ,  y  for  the  IP 

x*  y  xy 

be  compared  with  those  from  the  PB  in  Figs, 
mate  location  of  the  opened  crack  has  been 
orientation.  Moreover,  strain  contours  for 
15 f  should  be  ignored  since  small  strain 


The  strain  contours  are  generally  more  irregular  (especially  the 
shear  strain  ones)  for  the  IP  than  for  the  PB  and  that  was  expected  due 
to  the  presence  of  the  rigid  particles  in  the  IP.  However,  some 
irregularit ?es  are  seen  in  the  PB  strain  contours  as  well  which  may  be 
inherent  in  the  digitizing  process. 


Focusing  upon  a  comparison  of  the  e  contours. 


we  note  that  the 


intense  strain  zone  mead  of  the  blunted  crack  in  the  IP  (Fig.  5d)  tends 

to  increase  the  stma'n  ‘e/eis  surrounding  that  zone  above  those  in  the 

PB  (Fig.  51)  especbi'!/  i.i  the  direction  of  c^ack  extension.  A  similar 

conclusion  is  re.m  ^  I  in  ,.;m;  minq  Figs.  5e  and  5j  for  the  shear  strain 

contours  y  .  Tn-s  : .  ]<;•  ,ts  that  any  mathematical  model  which  is  used 
xy 

to  describe  the  simai"  f ieios  surrounding  the  near  tip  intense  strain 
zone  should  reflect  ’me  of  tnat  zone.  The  computer  plot  of  the 

contours  in  Fig.  3  i  for  the  PB  also  shows  an  intense  strain  field 


very  near  the  crack  tip.  nut  it  is  it  least  an  orde^  of  magnitude 
smaller  in  size  thm  in  ‘me  IP.  On  the  other  hard,  the  larger  crack 
opening  in  the  PB  rl  ‘me  rear  tip  blunting  in  the  IP  both  seem  to 
distribute  the  near  tic  sneer  strain  intensity  (Figs.  5e  and  5j)  over 
the  enlarged  tip.  •  ,  of  course,  van:shes  due  to  symmetry  at  the  crack 

*  i 


tip  on  the  crack  p '  --  •  . 
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Figure  6  shows  typical  crack  front  shapes  which  were  traced  from 
dye  markings  on  inert  propellant  and  pure  binder  respectively.  It  is 
clear  that  there  is  no  thumbnailing  in  the  inert  propellent  specimen. 
This  result  was  typical  of  all  tests  on  the  inert  propellant  and  «as 
somewhat  unexpected.  However,  if  one  focuses  upon  the  damage  2one  aneac 
of  the  crack,  one  may  conjecture  that  resin  rich  tubes  of  binder  which 
are  formed  here  (see  Figure  7)  may  not  be  strongly  transversely  con¬ 
nected,  thus  eliminating  the  transverse  constraint  normally  associated 
(as  a  precursor)  with  thumbnailing.  This  is  an  important  observation 
since,  if  true,  it  may  mean  that  no  plane  strain  fracture  toughness  will 
exist  for  this  material  in  this  environment. 

The  crack  fronts  displayed  by  the  pure  binder  were  even  more 
remarkable.  As  seen  in  Figure  6,  a  small  amount  of  reverse  thumbnailing 
occurred.  It  was  also  noticed  that  crack  growth  may  have  been  slightly 
retarded  by  the  grating  on  the  grating  side.  The  reason  for  the  inverse 
thumbnailing  is  not  clear  unless  the  surface  of  the  very  soft  binder 
material  could  have  been  damaged  in  making  the  initial  crack  and  this 
influenced  the  crack  front  shape  which  was  only  extended  by  about  2  mm. 

GRATING  SIDE 


DIRECTION  OF  GROWTH 


A 

PURE  BINOER  !83)  INERT  PROPELLANT  (PZ) 

Fig.  6  Crack  Border  Shapes  After  Some  Growth 


15.24  mm 


Fig.  7  Idealised  Material  Behavior  During  Crack  Opening  (large 

particles  spearate,  stretching  the  binder  into  tubes  between 
which  voids  develop)  (a)  No  Load,  (b)  Under  Load 

It  was  also  observed  that  severe  blunting  occurred  in  the  inert 
propellant  (as  noted  in  Ref.  1)  leading  to  the  use  of  a  modified 
algorithm  for  use  in  extracting  the  lowest  eigenvalue  \  from  the  data 
in  the  form:  u 

v-(v)0  ■ 

However,  in  the  pure  binder,  the  blunting  was  absent  and  (V)n  =  0. 
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(1) 


Figure  8  illustrates  the  use  of  the  equation: 


in( V-(V)Q !  =  jn  0^  +  xu2.n  r  (2) 


to  determine  \  .  Results  from  the  biaxial  tests  are  shown 

u 

in  Table  II. 

TABLE  2  -  Biaxial 

Test  Results 

Inert  Prop. 

Tests 

Test  PI 

Test  P2 

Linear  Zone 

(range) 

=  3. 5 -5.0  mm 

Linear  Zone 

(range) 

=  1. 5-5.0  mm 

Displ.  Pate 

=  2.54 

mm/min 

Displ.  Rate 

=  2.54  i 

mm/min 

a 

\ 

u 

2v^(mm) 

f% 

X 

u 

2vQ(mm) 

1.5% 

.629 

.102 

.76% 

.642 

0 

3.3% 

.744 

.476 

2.75% 

.674 

.20 

4.8% 

.532 

.143 

5 . 54% 

.582 

.16 

8.0% 

.648 

.344 

7.57% 

.652 

.26 

10.81% 

.792 

.13 

24.05% 

.712 

0.0 

vu(avg)  = 

0.638 

xu(avg)  =  0 

.676 

Test  P3 

Test  P4 

Linear  Zone 

(range) 

=  1.5-4. 5  mm 

Linear  Zone 

(range) 

=  1.5- 5.0  mm 

Displ.  Rate 

=  25.4  mm/min 

Displ.  Rate 

=  25.4  mm/min 

■:% 

X 

u 

2vQ(mm) 

t% 

X 

u 

2vQ(mm) 

.79% 

.592 

.044 

1.67% 

.643 

0 

1.58% 

.604 

.136 

5.64% 

.646 

.18 

3.16% 

.619 

.384 

6.43% 

.761 

0 

6.43% 

.570 

.247 

8.03% 

.677 

0 

10.38% 

.567 

.186 

13.64% 

.717 

.08 

12.74% 

.681 

.144 

13.54% 

.741 

.160 

xu(av9)  = 

0.625 

xu(avg)  =  0 

.689 
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TABLE  2  Biaxial  Tests  Results  Continued 


Pure  Binder  Tests 


Test  B1 


Test  B2 


Linear  Zone  (range)  =  2. 0-6.0  mm 
Oispl.  Rate  =  2.54  mm/min, 

2vq  =  0  mm 


Linear  Zone  (range)  =  1.5-5 
Displ.  Rate  =  2.54  mm/min, 
2vq  =  0  mm 


1.94% 

.727 

.37% 

.748 

3.49% 

.628 

1.93% 

.730 

5.06% 

.569 

3.47% 

.737 

6.60% 

.597 

12.87% 

.631 

9.74% 

.568 

17.36% 

.670 

17.35% 

.597 

\u(avg)  =  0.614 

xu(avg)  =  0.703 

Test  B3 

Test  B4 

Linear  Zone  (range)  =  1. 0-6.0  mm 

Linear  Zone  (range) 

=  1.5- 

Displ.  Rate  =  25 

.4  mm/min. 

Oispl .  Rate  =  25.4 

mm/min 

2vq  =  0  mm 

2vq  =  0  mm 

e% 

X 

e% 

X 

u 

u 

1.65% 

.807 

6.40% 

.749 

3.24% 

.642 

8.00% 

.722 

6.41% 

.569 

11.20% 

.660 

7.98% 

.573 

12.80% 

.668 

12.80% 

.526 

16.05% 

.556 

\u(avg)  =  0.612 

xu(avg)  =  .700 

Test  B5 

Test  86 

Linear  Zone  (range)  =  1. 5-5.0  mm 

Linear  Zone  (range) 

=  1.5- 

Displ.  Rate  =  2. 

54  mm/min, 

Displ.  Rate  -  2.54 

mm, 

2vg  =  0  mm 

2vq  =  0  mm 

€% 

\ 

e  % 

X 

u 

u 

4.68% 

.595 

6.24% 

.587 

7.80% 

.580 

12.52% 

.534 

xu(avg)  = 

0.574 

7.03% 

.641 

10.55% 

.608 

12.92% 

.556 

15.30% 

.566 

17.36% 

.602 

\u(avg)  =  0.595 


0  mm 


0  mm 


0  mm 
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Fig.  9  Crack  Tip  Profiles  for  Inert  Propellant  and  Pure  Binder 

The  current  results  for  the  inert  propellant  yield  values  of  \  at 

the  free  surface  which  are  similar  to  those  obtained  in  Ref.  1.  How¬ 
ever,  the  zone  dominated  by  these  values  is  further  from  the  crack  tip 
(0.5  to  1.5  mm  in  Ref.  1  vs.  1.5  to  5.0  mm)  herein.  The  similarity  in 
vu  values  would  be  expected  since  neither  set  of  crack  fronts  were 

thumbnailed.  The  shift  in  the  location  of  the  zone  where  x^  dominates 

is  also  expected  since  higher  global  strain  levels  (and  more  crack 
opening)  were  obtained  during  crack  growth  in  the  current  experiments 
than  in  the  measurements  prior  to  crack  growth  in  Ref.  1. 

Another  observation  is  that  the  zone  dominated  by  xu  was  as  far  as 

or  farther  from  the  crack  plane  in  the  pure  binder  tests  than  in  the 
current  inert  propellant  tests  even  though  Vg  =  0  in  the  former  tests. 

This  is  believed  due  to  the  larger  crack  opening  in  the  binder  tests  ac 
illustrated  in  Figure  9.  The  rather  large  scatter  in  the  \  values  from 

the  binder  tests  is  likely  partially  due  to  the  effect  of  the  reverse 
thumbnailing  which  caused  the  angie  of  intersection  between  the  free 
surface  and  the  crack  border  to  differ  significantly  from  90°. 

There  was  a  noticeable  difference  in  the  crack  growth  process 
observed  in  Ref.  1  and  that  observed  in  the  inert  propellant  in  the  cur¬ 
rent  program,  [n  the  tests  reported  in  Ref.  1,  the  crack  extended  by 
first  blunting  extensively  and  then  connecting  with  one  or  more  voids 
above  or  below  the  original  crack  plane,  producing  an  undulating  or 
meandering  path  about  the  original  crack  plane.  In  the  current  tests 
the  blunting,  although  still  present,  was  much  less  than  in  Ref.  1. 

This  appeared  to  focus  the  high  strain  zone  ahead  of  the  crack  closer  to 
the  original  crack  plane,  and  resulted  in  elimination  of  the  meandering 
effect.  This  is  described  pictorial ly  in  Figure  10.  The  severe 
blunting  in  Ref.  1  is  conjectured  to  be  due  primarily  to  resin  rich 
layers  of  material  on  the  free  surfaces  of  the  specimen  which  were 
absent  from  the  current  program. 

FROZEN  STRESS  TESTS 

Since  thumbnailing  did  not  occur  along  the  crack  front  in  the 
biaxial  tests  on  the  inert  propellant,  the  frozen  stress  data  for  esti¬ 
mating  the  variation  of  \  with  thickness  were  taken  from  four  point 

bend  tests  on  straight  front  cracks  in  a  stress  freezing  epoxy  material. 
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M  F  =48 

The  specimen  dimensions  are  given  in  Figure  11.  A  companion  test  series 
employing  naturally  grown  cracks  which  did  thumbnail  was  also  conducted. 


P  <  0  050  mm 


W  =  25  70  mm 
d  =  6.35  mm 
a  =  12.82  mm 
B  =  13  33  mm 
e  =  12.7  mm 


I  =  50  8  mm 
L  =  762  mm 
S  =  279  4  mm 
P  =  2  32  N 


Fig.  11  Four  Point  Bend  Specimen  Dimensions 

All  such  specimens  were  prepared  from  stress  free  di-phase  trans¬ 
parent  epoxy  plates.  Fairly  dense  moire  gratings  (600  and  1200  lines/mm 

1.75x10  mm  thick)  were  cemented  to  the  free  surface  in  the  neighbor¬ 
hood  of  the  point  of  intersection  of  the  crack  front  with  the  free  sur¬ 
face  and  the  specimens  were  then  placed  in  a  cam-controlled  oven  and 
heated  to  critical  temperature  (=  225 ’F  or  107"C).  After  a  thermal 
soak,  four  point  loads  were  applied,  and  the  natural  cracks  were  grown 
to  desired  size.  After  cooling  under  load,  specimens  were  removed  from 
the  oven  and  the  deformed  gratings  at  the  free  surface  were  viewed 
through  a  master  virtual  grating  which  was  made  by  splitting  a  laser 
beam  and  reflecting  one  part  back  into  the  path  of  the  other  part  to 
form  walls  of  constructive  and  destructive  interference.  This  virtual 
grating  was  tuned  to  twice  the  frequency  of  the  active  grating  on  the 
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specimen  surface  with  a  diffraction  grating.  A  typical  moire  pattern 
for  the  displacement  component  normal  to  the  tip  of  a  thumbnailed  crack 
is  shown  in  Figure  12.  After  obtaining  vu  from  this  pattern  using  Eq.  ? 

(with  Vq  =  0),  slices  were  removed  mutually  orthogonal  to  the  crack 

front  and  the  crack  plane,  and  analyzed  photoelastically.  Using  the 

ir 

variable  x  algorithm  described  in  Appendix  B  (Eq.  B-7),  values  for  x^ 

were  obtained  through  the  thickness  of  the  specimens.  Figure  13 
compares  the  distributions  of  x^  for  the  straight  front  and  thumbnailed 

cracks.  The  x  distributions  from  two  thumbnailed  cracks  whose  crack 

a 

fronts  are  also  shown  in  Fig.  13  are  superimposed  on  the  distribution 
for  two  straight  front  cracks,  the  latter  results  taken  from  Ref.  1. 
Considering  the  inevitable  data  scatter  present  in  such  data,  one 
concludes  that  any  differences  between  the  two  data  sets  appear  to  be  of 
the  order  of  the  data  scatter  itself.  Part  of  the  scatter  may  be  due  to 
lack  of  symmetry  of  the  natural  cracks  but  no  adjustments  in  the  data 
were  made  to  account  for  this  effect. 


M.F.  =  18.8 


Fig.  12  Moire  Pattern  for  Near  Tip  Displacement  Field  (V)  at  Free 
Surface  of1  Inert  Propellant  Containing  a  Natural  Crack 


Figure  14  shows  a  comparison  between  the  normalized  distribu¬ 
tions  for  the  straight  front  and  natural  cracks  computed  from  Eq.  (B-4) 
in  Appendix  B.  The  "hump"  in  the  distribution  for  the  thumbnailed 
cracks  occurs  near  where  the  curvature  of  the  thumbnailed  cracks 
reverses  and  may  be  related  to  that  event  (marked  with  tick  marks  in 
Fig.  13).  Moreover,  growth  of  the  thumbnailed  cracks  was  not  symmetri¬ 
cal.  Slices  were  taken  both  parallel  to  the  free  surface  and  normal  to 
the  crack  front.  While  the  values  of  Kj  were  always  higher  for  slices 

taken  normal  to  the  crack  front  (as  expected),  there  was  less  data 
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Fig.  13  x  Distributions  for  Straight  Front  and  Thumbnailed  Cracks 


scatter  in  the  slices  taken  parallel  to  the  free  surface.  Both  slicing 
methods,  however,  revealed  the  hump  in  the  distribution.  Near  the 
hump,  scatter  was  of  the  order  of  i  6%.  Although  part  of  this  scatter 
could  be  due  to  optical  and  material  noise,  it  is  more  likely  the  result 
of  crack  front  curvature  which  was  not  symmetrical  with  respect  to  the 
center  of  the  thickness. 


SUMMARY 

Optical  methods  which  exert  negligible  influence  on  near  tip 
behavior  of  both  inert  propellant  and  pure  binder  in  through  cracked 
bodies  were  used  to  measure  the  displacement  fields  on  the  free  surface 
of  biaxial  test  specimens  in  the  neighborhood  of  the  crack  front-free 
surface  intersection  point  under  Mode  I  loading  at  head  rates  differing 
by  an  order  of  magnitude.  For  this  purpose,  grating  methods  were 
employed.  Data  collected  were  digitized  and  surface  in-plane  displace¬ 
ment  and  strain  fields  were  plotted.  The  value  of  the  dominant  eigen¬ 
value  at  the  free  surface  was  extracted  from  the  displacement  field 

measurements. 
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Fig.  14  Stress  Intensity  Factor  Distributions  in  Four  Point  Bend 
Specimens  with  Straight  and  Thumbnailed  Crack  Fronts 


Using  the  results  from  high  density  moire  interferometry  and  frozen 
stress  analysis  of  four  point  bend  specimens  of  di-phase  photoelastic 
materials,  the  distribution  of  \  through  the  thickness  was  inferred  for 

both  straight  front  and  "thumbnai led"  cracks. 

The  studies  described  above  together  with  those  of  Ref.  1  lead  to 
the  following  observations  and  conclusions  regarding  near  tip  behavior 
of  inert  propellant  in  cracked  biaxial  specimens  under  Mode  I  loading. 

OBSERVATIONS  AND  CONCLUSIONS 

The  crack  extension  process  in  inert  propellant  consists  of  two 
phases;  firstly,  crack  opening  and  blunting  during  which  severe  stretch¬ 
ing  of  the  binder  material  ahead  of  the  crack  produces  a  binder-rich 
craze-like  zone  of  binder  tubes  containing  voids  (Fig.  7)  ahead  of  the 
crack.  The  blunting  phase  is  due  to  the  presence  of  the  rigid  particles 
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of  inert  propellant  and  did  not  occur  in  tests  on  pure  binder.  The 
direction  of  crack  extension  may  oscillate  about  the  plane  of  the  main 
crack  if  excess  binder  is  present  (Fig.  10). 

When  the  cracks  grow  in  inert  propellant,  the  crack  front  remains 
straight,  indicating  an  absence  of  transverse  constraint  in  the  craze 
zone  ahead  of  the  crack.  This  is  an  important  observation,  since  it 
means  that  a  plane  strain  fracture  toughness  may  not  exist  for  this 
material  at  room  temperature. 

Values  of  the  lowest  dominant  eigenvalue  can  be  accurately 

measured  at  the  free  surface  of  the  inert  propellant  provided  the 
influence  of  blunting  is  properly  accounted  for  in  the  algorithm  for  the 
displacement  field  (Eq.  2). 

The  four  point  bend  test  results  appear  to  yield  a  good  description 
of  the  variation  of  \  from  its  free  surface  value  back  to  the  value  of 

one  half  as  per  LEFM  (Fig.  13). 

Changing  the  crosshead  speed  from  2.5  to  25.4  mm/min.  altered  the 
displacement  and  strain  fields  significantly  but  not  the  global  contour 
patterns  or  values  of  \  . 

When  compared  to  results  from  cracked  test  specimens  of  pure 
binder,  the  main  influence  of  the  rigid  particles  in  IP  on  near  tip 
behavior  is  to  enlarge  an  intense  zone  of  normal  strain  normal  to  the 
crack  plane  to  a  size  at  least  an  order  of  magnitude  greater  than  its 
size  in  cracked  specimens  of  pure  binder.  This,  in  turn,  enlarges  the 
surrounding  strain  fields  (Fig.  5).  The  rigid  particles  also  served  to 
strengthen  and  stiffen  the  material  globally  (Fig.  3). 
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APPENDIX  A 


DIGITIZATION  PROCESS 

Points  located  at  the  intersection  of  the  grating  lines  (such  as 
those  shown  in  Fig.  4)  are  selected  close  together  in  regions  of 
expected  high  gradients  and  further  apart  away  from  these  regions. 

These  points  in  groups  of  four  form  quadri laterals  (initially  rectan¬ 
gular)  and  values  are  read  at  each  of  four  points  and  averaged  before 
and  after  loading.  The  four  point  average  is  located  at  the  centroid  of 
each  quadri lateral  and  the  difference  between  the  no-load  and  loaded 
values  becomes  the  digitized  displacement  at  that  point.  An  array  of 
these  Tatter  points  before  and  after  loading  are  shown  in  Figs.  A- 1  and 
A-2  along  with  the  correspond ing  displacement  contours  plotted  by  com¬ 
puter.  The  computer  program  plots  contours  representing  the  displace¬ 
ments  of  the  points  from  their  initial  positions  (Fig.  A-l)  to  their 
final  positions  (Fig.  A-2). 
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Fig.  A-2  Points  for  Digitization  After  Loading 
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APPENDIX  B 


LEFM  FROZEN  STRESS  ALGORITHM  -  TWO  PARAMETER  APPROACH 

By  choosing  a  data  zone  sufficiently  close  to  the  crack  tip  that  a 
Taylor  Series  Expansion  of  the  non-singular  stresses  can  be  truncated  to 
the  leading  terms,  one  may  write  for  a  straight  front  crack  [B-l| 


- -  tto  .  iJ  =  1.2 


(B-l) 


where  a.  .  are  the  stress  components  in  the  plane  normal  to  the 

crack  border,  is  the  Mode  I  Stress  Intensity  Factor  and  r,  9  are  the 

polar  coordinates  with  origin  at  the  crack  tip  (see  Fig.  B-I).  Along 

0  2  2 

9  =  n/2,  one  has,  for  case  where  (0  )  is  small  relative  to  8-t  (see 

0  m 
Ref.  I B- 1 ]  (where  o  is  proportional  to  the  contribution  of  the 

nonsingular  stress  to  the  maximum  shear  stress,  imax  in  the  x^x^  plane) 


Now  defining  an  "apparent"  SIF, 


KflP  *  'max(8"r)1/Z  <B-3> 

-  1/2  - 

and  normalizing  with  respect  to  o(na)  where  0  represents  the  remote 
stress  and  a  the  crack  length,  one  has 
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a(  Tra) 


Ki  0  .  n 

= _ 1 _  4-  2_  fl)1/2 

- .  1  n  ~  U  J 


(B— 4) 


which  suggests  an  elastic  linear  2one  (ELZ)  in  a  plot  of 


r  1/2  u 

(-)  '  with  a  slope  of  —  .  Experience  shows  this  zone  to  lie  between 

d 

(r/a)^2  values  of  approximately  0.2  to  0.4  (or  above)  in  two  dimen¬ 
sional  problems.  By  extracting  optical  data  from  this  zone  and  extra¬ 
polating  across  a  near  tip  nonlinear  zone,  an  accurate  estimate  of 
-  1/2 

K^/rt^a)  can  be  obtained.  This  is  illustrated  in  Figure  |B-2|  for 

two  slices  at  different  locations  in  a  four  point  bend  specimen.  This 
approach  may  be  directly  extended  to  3D  problems  involving  curved  crack 
fronts  by  replacing  i,j  values  of  1,2  by  n,z  where  the  nz  plane  is  nor¬ 
mal  to  the  crack  border  at  each  point  along  the  flaw  border  and  the 
n,z,t  coordinate  system  is  a  local  cartesian  system  which  follows  the 
crack  border  (Figure  8-1).  In  three  dimensional  problems  the  outer 

boundary  of  the  data  zone  is  usually  restricted  to  (r/a)^  =  0.4  or 
less. 

Comparison  of  Photoelastic  Doto  for  Interior 
and  Exterior  Slices 

H — DEPTH — -j 


Kth  from  A  STM  E399  (ZD) 


Fig.  B-2  Determination  of  Stress  Intensity  Factors  from  Frozen  Slice 
Data  in  Four  Point  Bend  Specimen 
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VARIABLE  Xo  ALGORITHM 

Following  the  above  approach  and  the  work  of  Benthem  Ref.  [B-2i, 
one  may  formulate  an  expression  for  the  near  tip  o.^  in  the  form: 

^  i  j  ^  X  3  *  6  ^  0  r 


OJ  =  n,z) 


(B-5) 


where  a9^  =  for  i,j  =  n  and  zero  for  all  other  values  of  i,j.  From 

Eqs.  (B-5)  one  may  develop  an  algorithm  from  which  values  of  x^  may  be 

extracted  but  the  procedure  is  not  as  straightforward  due  to  the 
presence  of  the  second  term  in  Eq.  B-5  and  the  highly  three  dimensional 
nature  of  the  near  tip  stress  state  at  the  free  surface.  If  one 
evaluates  r  from  Eq.  (B-5)  along  8  =  tt/2  as  done  in  LEFM,  one 
obtains:  max 


7  *  r  s1n<'0  ♦ ')  ?  ■ 


/2i  r 


«x  >, 


/2-t  r 


(B-6) 


where  K  may  be  designated  a  stress  eigenfactor.  Defining 

a°  0 

sin(x_j  +  1)  ^  ,  one  obtains 


in^max  -  V  =  ln  *r-=?:  “  Vln  r> 

/d  TT 


(B-7) 
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*  tq  is  first  computed  from  LEFM  at  all  points  (i.e.  using  x^  =  1/2) 

except  at  the  free  surface  where  we  set  =  0.  Then,  at  the  onset 
of  decreasing  Tg  below  midthickness  values,  we  decrease  linearly 
from  that  point  to  zero  at  the  free  surface.  This  approach 
adequately  accounts  for  the  effect  of  the  nonsingular  stress. 
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